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We calculate the uj — > it 7* and <f> — > 7r 7* electromagnetic transition form factors based on 
dispersion theory, relying solely on a previous dispersive analysis of the corresponding three-pion 
decays and the pion vector form factor. We compare our findings to recent measurements of the 
uj Tr°jj, + n~ decay spectrum by the NA60 collaboration, and strongly encourage experimental 
investigation of the Okubo-Zweig-Iizuka-forbidden <f> —¥ ir°£ + £~ decays in order to understand the 
strong deviations from vector-meson dominance found in these transition form factors. 
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INTRODUCTION 



In recent years, there has been intense renewed interest 
in light-meson transition form factors due to their poten- 
tial role in the theoretical determinations of the anoma- 
lous magnetic moment of the muon (see Ref. |l| for a re- 
view) . With more and more exclusive channels contribut- 
ing to the hadronic vacuum polarization measured ex- 
perimentally with unprecedented precision, it is believed 
that the hadronic contribution to light-by-light scattering 
may soon constitute the dominant uncertainty 2] . While 
a full determination of the light-by- light scattering tensor 
remains a formidable task, a combination of experimen- 
tal data and theoretical analyses may help to constrain 
one of the most important contributions (and one of the 
few that are model- independently accessible), namely the 
pseudoscalar (P = 7r°, 77, rf) pole terms. Their strength 
is determined e.g. for the ir° pole contribution by the de- 
cay 7T° — > 7*7*, given in terms of the doubly- virtual form 
factor F T o 7 » 7 . (M% , <?|) ( see e -§- R- ei - PJ f° r precise 
definitions), where q^ 2 denote the two photon virtual- 
ities. As these doubly- virtual form factors, that are to 
be measured in the rare decays P — > £ + £~£' + £'~ (with 
branching ratios of the order of 10 -5 ), are difficult to de- 
termine precisely in experiment, it is useful to note that 
they are intimately linked (for specific values of one of the 
photon virtualities) to vector-meson conversion decays: 
e.g., the form factor F 7r o 7 , J ,(M^ ,q 2 ,M 2 ) determines 
the dilepton spectrum in w — > ir°£ + £~, the form factor 



F 
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can be measured in — » rj£ + £ etc. 



The interactions of hadrons with (real and virtual) 
photons are often thought to be described at least to 
good approximation in the picture of vector-meson dom- 
inance (VMD): the q 2 -dependence of the form factors 
above should largely be given by the propagator of a light 
intermediate vector meson (ft, uj, </>), see e.g. Refs. 
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What is interesting about the vector-meson conversion 
decays is that they show a very clear deviation from such 
a simple VMD picture, as has been established in the de- 
cay uj — > 7r°/i + /i~ 0-01, and also in <j> — > rje + e~ [ToL [TTI| . 

In this article, we will analyze two such vector-meson 
transition form factors with the method of dispersion 
relations, concentrating on uj — > 7r°7* and <fi — > 7r°7* 
(the latter being rarer due to the implied violation of 
the Okubo-Zweig-Iizuka rule; see Ref. [l2[ for a recent 
theoretical work). One specific theoretical advantage of 
vector-meson conversion decays, as opposed to the pseu- 
doscalar Dalitz decays, is that the isospin of the virtual 
photon is fixed (in the approximation that isospin is con- 
served). In the cases at hand, it needs to be an isovector 
photon, hence the lowest-lying intermediate states to con- 
tribute in a dispersion relation are 2ir, 4tt etc. only, and 
experience with pion-pion P-wave interactions suggests 
that the 2tt intermediate state will already saturate the 
dispersion relation to a large degree. As we will demon- 
strate below, a dispersive reconstruction of the 2ix contri- 
bution requires two amplitudes as input: the correspond- 
ing V — > 7T + 7r~7r° decay amplitude (in the appropriate 
partial wave), and the pion (electromagnetic) vector form 
factor. 

An analysis of the uj — » 7r°7* transition form factor us- 
ing dispersion theory has already been performed decades 
ago [13[ (although phenomenologically the focus of that 
work lay more on the e + e~ — > ujit production cross sec- 
tion) . The reasons to take up this subject again are mani- 
fold: we now have much more accurate experimental as 
well as theoretical input at our disposal, both for the 
pion vector form factor and the required pion-pion phase 
shifts; furthermore, we have recently performed a disper- 
sive analysis of the three-pion decays of both uj and <f> that 
treats final-state interactions between all three pions rig- 
orously [lij , which can now serve as the consistent input 
to the investigation of the transition form factors. All 
these ingredients of the analysis will be reviewed below. 

As a final introductory remark, we note that in this 
article, we will confine ourselves to an analysis of the 
transition form factors in the kinematical region acces- 
sible in the corresponding vector-meson decays. We are 
aiming for a precision analysis and hence refrain from 
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analyzing also the processes e + e~ — > tt°uj |15l - [l7| and 
e+e~ — > 7r°0 [18], as we expect these to be signifi- 
cantly more dependent on information from the excited- 
resonance region. 

The outline of this article is as follows. We intro- 
duce the necessary definitions concerning kinematics and 
partial- wave decomposition in SecHT] We discuss the dis- 
persion relation for the transition form factors in Sec. Mil 
including the two main elements required as input: the 
pion vector form factor and the V — > ir + ir~ir° partial- 
wave amplitude. Numerical results for form factors, 
decay spectra, and branching ratios are presented in 
Sec.|lVl before we summarize in Sec. El Technical details 
on a representation of the pion vector form factor includ- 
ing higher resonances are relegated to an Appendix. 



II. KINEMATICS AND PARTIAL- WAVE 
DECOMPOSITION 

We consider the decays of the lightest isoscalar vector 
mesons into a ir° and a dilepton pair, 



V(p v ) -> 7r°(p )£ + {p e+ )e-( Pe -) , V = u/4>, I = e/fi. 



The V — > tt°£ + 1 amplitude is given as [7j| 



(1) 



My-* = ie e^apn^poq 



u s (pe-hv S '(pt+) 



(2) 

where q = pi+ + Pe-, s = (pv — Po) 2 , is the polar- 
ization vector of the vector meson, and fv-n a ( s ) is the 
electromagnetic transition form factor of the vector me- 
son. We will also discuss the corresponding normalized 
form factor, 



fvirO Q) 



(3) 



The differential decay rate in terms of this amplitude can 
be written according to 



dT 



v^-KH+t- 2a Vi , 2m2 i\ <?" 0)<?vv° ( s ) If ( M2 
-Js = + — ™ I/ ^° (S)I ' 



M y 



(4) 

with the fine-structure constant a = e 2 /4ir and masses of 
vector meson, neutral pion, and leptons denoted by My, 
M n o, and mi, respectively. The center-of-mass momenta 
are given by 



X(M 2 A ,M%,s) 



4s 



(5) 



where A(x, y,z) — x 2 + y 2 + z 2 — 1(xy + yz + xz) is 
the Kallen function (and with the slight notational abuse 
Mf. = rag implied). Radiative corrections to Eq. Q have 
been calculated in Ref. [20|: they require a careful selec- 
tion of kinematic cuts on the additional soft-photon radi- 
ation for the e + e~ final state, and are small everywhere 
except near threshold for fi + fi~, where the Coulomb pole 




FIG. 1: Diagrammatic representation of the discontinuity of 
the V — > tt°£ + £~ transition form factor. The gray circle de- 
notes the V — > 37r amplitude, whereas the white circle repre- 
sents the pion vector form factor. 



is significant. The above relation for the V — > tt°1 + 1~ 
spectrum is completely determined by fy^o (s) aside from 
a kinematical factor that is determined by the photon 
propagator and phase space. Note finally that the corre- 
sponding real-photon total decay rate is given by 
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V— >7T°7 



a(M v -M 2 ) 3 
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In establishing a dispersion relation for the V 
transition form factors, the corresponding three-pion de- 
cays V(pv) — > tt + {p+)tt~ (p-)vr (p ) play a central role. 
We define s — (pv — Po) 2 , t — {pv ~ P+) 2 , and u — 
(p v ~ P-) 2 with 3s = s + t + u = M v + 3M 2 . We 
note that due to technical reasons, the decay V — » 3ir 
is treated in the isospin limit, with M v o = M^± = M v . 
The amplitude is given as 



M 3 * 



^p v +P 1p P T(s,t,; 



(7) 



Neglecting the discontinuities of F and higher partial 
waves, 1 we can decompose J-(s, t, u) in terms of functions 
of a single variable with only a right-hand cut flij ]. 

T(s, t, u) = T(s) + F(t) + F(u) . (8) 

For our analysis we will require the / = 1 partial-wave 
projection of T(s, t, it), which is given by 



his) = 



3 f 1 

— J dz(l — z 2 ) J-(s,t,u) 



(9) 



where z = (t — u)/(4g ffw (s)qy 7r o(s)). Note that the angu- 
lar integration in Eq. © is a highly nontrivial issue due 
to the fact that the cuts in the variables t and u need 
to be avoided; compare the discussion in Appendix B of 
Ref. [JJl. 



III. 



DISPERSION RELATION FOR THE 
TRANSITION FORM FACTOR 



To set up the dispersion relation for the transition form 
factor, we calculate the two-pion discontinuity of the di- 



1 A simplified model for additional F-wave contributions was stud- 
ied in Ref. Il4l and found to yield entirely negligible corrections. 
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agram shown in Fig. [TJ It is given as [R 



disc f Vv o (s) 



lq ^ is) -FY*(s)h( S )9( S - 



4Mi 



(10) 



where Fjf(s) is the pion vector form factor. Corrections 
to Eq. (fTU|) stem from heavier intermediate states of the 
appropriate quantum numbers (isospin 1 P-wave states): 
Att, K K, . . . , which are expected to be suppressed signif- 
icantly due to phase space and their higher masses. We 
therefore neglect these contributions in our analysis and 
resort to elastic tttt final states only. Given the standard 
assumptions on the asymptotic high-energy behavior of 
the pion form factor, (s) ~ 1/s (modulo logarithms), 
and the V -> 3?r partial wave, / x (s) ~ 1/s [lj], Eq. (fTTTj) 
allows for an unsubtracted dispersion relation [13]. As 
our analysis, however, is confined to two-pion interme- 
diate states and neglects any higher contributions, we 
decide to employ a once-subtracted solution of Eq. (fit)]) 
instead, 



f v „o(s) = /y T o(0) + 



12tt 2 



ds 



,q^{s')Fy*{ S ')h{ S ') 



4MI 



s /3 / 2 (s'-s) 



in order to suppress inelastic contributions. For the pre- 
dictions of the s-dependence of the form factor, we fix the 
subtraction constant fvn (0) to reproduce the V — > tt°j 
partial width according to Eq. ([6]) . Assuming the validity 
of an unsubtracted dispersion relation, the subtraction 
constant and therefore the V — > tt°j partial width can 
be calculated by means of a sum rule. 



iWo(o) 



i 



12tt 2 



ds 



/<77T7r( S ') 77.V* 



IM'l 



s /3/2 



FV*(s')h(s') , (12) 



which is expected to show a more problematic conver- 
gence behavior than the form-factor dispersion relation. 
We will quote results for this sum rule in Sec. IIVI in or- 
der to quantify the potential role of heavier intermedi- 
ate states in the transition form factor. Nevertheless, 
Eq. (TT2"j) is a remarkable result: as we will briefly reiter- 
ate below, in the elastic approximation the pion vector 
form factor is entirely given in terms of the tttt P-wave 
phase shift £>f(s), which also determines the V — > 3tt par- 
tial wave up to a single subtraction constant that can be 
written as an overall normalization [l4| . This means that 
the ratio of branching ratios B(V — > tt®^)/B{V — > 37r) is 
entirely determined by 5\{s), up to inelastic corrections. 
This result is reminiscent of the relation between these 
two decay modes utilized in Refs. [U [22[ as the leading 
order of a Lagrangian framework for vector mesons. 

In the following we will briefly discuss the two ingre- 
dients to the dispersion integral, the pion vector form 
factor F^(s) and the V — ¥ 37r partial-wave amplitude 
Ms). 
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FIG. 2: Pion vector form factor fit using a phase shift incorpo- 
rating elastic p', p" resonances (solid line) and solutions of the 
Roy equations of Refs. 24, 2{| (dashed and dotted lines), in 
comparison to the experimental data of Ref. [2(| . For details, 
see main text. 



A. Pion vector form factor 

In the elastic approximation, the pion vector form fac- 
tor fulfills the unitarity relation 

disc-Fj(s) = 2iF^(s)9(s - 4M^)smd(s)e- iS{s) , (13) 

where S(s) = 8{(s) is the tttt P-wave phase shift. The 
solution to Eq. (|13p is given by the Omnes function, 



F?(a) = Sl(8) 



cxp< - 



S(s') 



ds' 

AMI s '\ s> - s ) 



(14) 



normalized to O(0) = 1. The omission of a polynomial in 
s multiplying the Omnes function relies on the absence 
of zeros in the form factor, see Ref. [23]. In a precision 
analysis of the form factor extending beyond 1 GeV, one 
has to account for the onset of inelasticities (dominantly 
Att intermediate states), and, as far as data extracted 
from e + e" — > tt + tt~ is concerned, p-uj mixing. As we do 
not have a consistent treatment of inelasticity effects in 
the V — > 3tt partial wave /i (s) at our disposal (let alone 
isospin breaking), we refrain from doing so. 

We use the following approach to estimate uncertain- 
ties generated by the input for the parameterization of 
the phase shift. On the one hand we will use param- 
eterizations derived from two different solutions of the 
pion-pion Roy equations [24], [25| , which are valid roughly 
up to 1.3 GeV. The pion form factor is known to excel- 
lent precision well beyond that energy (see Refs. [26M28| 
for just the most recent experimental results), indicating 
in particular contributions from the excited resonances 
p'(1450) and p"(1700). To incorporate these higher res- 
onance states we use the phenomenological form factor 
suggested in Ref. (29j (and briefly summarized in Ap- 
pendix which we fit to the experimental data of 
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Ref. [26j, extract the corresponding phase, and match 
it smoothly to the phase-shift solution of Ref. [24J below 
1 GeV. The aforementioned procedure treats the higher 
resonances p' and p" as purely elastic, which they clearly 
are not (compare the more sophisticated form factor rep- 
resentation of Ref. 30]); we merely use the phase thus 
obtained as an indicator for uncertainties generated in 
the energy range between roughly 1.3 GeV and 1.9 GeV. 

As the Omnes representation requires the tttt P-wave 
phase shift up to infinity, we have to make assumptions 
about its asymptotic behavior. We choose to smoothly 
guide S(s) to ir, so that we guarantee the correct asymp- 
totic behavior of (s) —> s _1 for s —> oo. The point 
beyond which the asymptotic behavior sets in is chosen to 
be As = 1.3 GeV for the Roy-equation analyses [HHHi, 
and A,5 = 1.9 GeV for the phase derived from the form 
factor. The resulting form factors corresponding to the 
different phases are shown in Fig. [5J compared to the 
data of Ref. 



B. V 



3tt partial-wave amplitude 



We only very briefly summarize the basics of the dis- 
persive analysis of the V — > 3tt decay amplitudes and 
refer to Ref. [l4| for further details. The partial-wave 
amplitude /i(s) fulfills the unitarity relation 

disc /i (s) = disc J"(s) (15) 
= 2t (.F(s) + F(s))9(s - AMI) sin6{s)e~ u{s) , 

where the inhomogeneity F(s) is given by angular aver- 
ages over T according to 



(16) 



(z n f)(s) 
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3s - s 
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The function -F(s) contains the left-hand-cut contribu- 
tions due to crossed-channel singularities. In the case 
at hand, the left-hand cut overlaps with the right-hand 
one, as for My > 3M W , s, t, and u can be simultaneously 
larger than 4M^, which they are inside the physical de- 
cay region. In this sense, there is no meaningful way 
from the point of view of dispersion theory to neglect the 
left-hand cut here. The solution of Eq. ([15")) is given by 
a once-subtracted dispersion relation, 



T(s) = Q(s)\ a + 



da' sin S(s')f(s') 
tt J iM , ~V \n( s >)\( s > - S ) 



(17) 



where the subtraction constant a serves as an overall nor- 
malization and is adjusted to reproduce the V — ¥ Sir par- 
tial width. 

As Ref. [HI does not explicitly show the partial-wave 
amplitude /i(s) = J-(s) + T(s) that plays a central role 
in the present investigation, we display its modulus and 




FIG. 3: Modulus (upper panel) and phase (lower panel) of the 
P partial wave /i (s), based on Ref. [1J| both for uj — » 37r (solid 
line) and — > 37r (dotted line), in comparison to the Omnes 
function (dashed line). We refrain from devising error bands 
and fix the input for the phase according to Ref. [24[ and the 
integral cutoff in Eq. (|17|l to A = 2.5 GeV. The normalization 
constant a is set to 1. 



phase both for u — > 3n and & — > 37r, as derived from the 
numerical results in Ref. [lj], compared to the Omnes 
function (whose phase of course is just 8(s)) in Fig. [3] 
We note that the partial waves bear very little similarity 
to the Omnes function: there is a strong enhancement in 
the threshold region below the p resonance, a large part 
of which can be thought of as the partial- wave- pro jected 
t- and u-channel p exchanges in a VMD picture [13fl . Fur- 
thermore, we note that the phase of /i(s) also does not 
follow S(s): Watson's theorem does not hold due to three- 
pion-cut effects, see Fig. [H which in particular allow for 
a nonvanishing imaginary part of fi(s) already at tttt 
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FIG. 4: Two-loop diagram contributing to the V — ¥ 3ir decay 
amplitude, which has a singular discontinuity at the pseudo- 
threshold s = (Mv — Mtt) 2 and leads to a nonvanishing phase/ 
imaginary part of the corresponding partial wave /i(s) at 
threshold s = 4M 2 . 



threshold. 2 Finally, /i(s) shows singular behavior at the 
pseudothreshold s — (My — M^) 2 , stemming from the 
diagrammatic topology shown in Fig. HJ as discussed in 
Refs. dHHH, these singularities in the discontinuity do 
not translate into singular behavior of the form factor 
itself when evaluated at the upper rim of the unitarity 
cut. In particular, the irregular phase in the vicinity of 
the pseudothreshold is an artifact as a consequence of the 
different divergences of real and imaginary parts from be- 
low and above, and has no physical significance. 

Returning to the dispersive representation of ^(s), it 
was shown in Ref. [lj] that by oversubtracting the inte- 
gral equation (flTl) for 4> — > 37r, 



bs + — 

7T 



4M2 \n( s >)\(s> - s) 



(18) 

and adjusting the additional subtraction constant b, we 
were able to achieve a perfect fit of the </>—>• 3 ir Dalitz 
plot H. Since the fitted value for b differs from a 
sum rule, as suggested by demanding the representa- 
tions (fTT)) and (JTHJ) to be equal, Eq. ([T%1) does not sat- 
isfy the high-energy behavior for the partial- wave ampli- 
tude /i(s), which therefore tends asymptotically towards 
a constant instead of s _1 ; consequently, the integral (fT2")l 
does not converge, and we will not evaluate the sum rule 
for /^(O) for the twice-subtracted solution of fi{s). 

We will in general stabilize the high-energy behavior 
of our dispersion integrals by manually leading T(s) to 
A 2 J r (A 2 )/s beyond a certain cutoff A 2 . There is no ob- 
vious prescription as to when exactly the amplitude or 
the form factor should show this asymptotic behavior; we 
choose the point up to where we have adjusted our form 
factor representation to data, that is A = 1.8 GeV, and 
incorporate a variation of the cutoff up to A = 2.5 GeV 



2 Note that this complication does not occur in Ref. [13j due to 
the approximation of the left-hand-cut contributions by p poles 
only, neglecting the effects of the two-pion cut starting already 
at t, u = 4M 2 .. As only the transition form factor of the lj 
and not the one of the tf> is considered in Ref. |13( . the p poles 
in t-/M-channel lie outside the integration range that affects the 
partial-wave projection. 



in our error considerations. This prescription assures 
that we have a precision representation for the ampli- 
tude in the low-energy regime as well as the correct high- 
energy behavior. By varying the cutoff we assure that 
the intermediate-energy regime is sufficiently suppressed 
so as not to taint our numerical results, which we present 
in the following section. 



IV. NUMERICAL RESULTS 

For the numerical evaluation, we use the different pa- 
rameterizations of the phase shift described in Sec. IHI Al 
the same parameterization is always used consistently for 
both pion form factor and V — > 3tt partial wave. We 
vary the cutoff of the dispersion integrals in Eqs. (|lll) . 
(f!2l) . and ([l~7|) (beyond which the assumed asymptotic 
behavior is enforced by hand) between A = 1.8 GeV and 
A = 2.5 GeV as detailed above. We note that it does 
not make sense to vary the cutoff of the aforementioned 
integrals individually: the uncertainties in our treatment 
are related to our lack of knowledge concerning final- 
state interactions in the intermediate-energy range, and 
thus apply equally to all considered dispersion integrals. 
The subtraction constants of the V — ¥ 3n amplitudes 
arc fixed by the total widths = 8.49 ± 0.08 MeV and 
r = 4.26 ±0.04 MeV together with the V -> 3tt branch- 
ing ratios £ exp (w -> 3tt) = 0.892±0.007, B cxp {cj) -> 3tt) = 
0.153 ± 0.003, the uncertainties of which we will always 
neglect in the following. It turns out that for all of our 
evaluations of the once- subtracted dispersion relation in 
Eq. (jlip . a variation of the phase between the parameter- 
ization of Ref. [25] and the one derived from the form fac- 
tor spectrum along with an integral cutoff of A = 1.8 GeV 
gives rise to an enveloping uncertainty band. 

For the following V —> 7r°7 branching ratios deter- 
mined from Eq. ([6]) and the unsubtracted dispersion re- 
lation (fT2"j) , the parameterization of Ref. [24J and the one 
from the form factor spectrum together with an integral 
cutoff of A = 2.5 GeV give rise to limiting values. We 
find 



B{lj -> tt°7) = (7.48. 
S(0^7r° 7 ) = (1.28. 



.7.75) x 10" 2 
.1.37) x 10~ 3 



(19) 



which is to be checked against the experimental aver- 
ages B cxp (lj -> tt°7) = (8.28 ± 0.28) x 10~ 2 , £ cxp (0 -> 
7T°7) = (1.27 ± 0.06) x 10" 3 HI. We observe that the 
(f> — > 7r°7 partial width compares favorably to experi- 
ment, whereas the result for oj 7r°7 turns out to be 
slightly too low; even then, the 27r intermediate state 
seems to saturate more than 90% of the sum rule for this 
partial width. We note that the most precise individual 
measurement of B(oj — > 7r°7) actually determines the ra- 
tio of branching ratios B exp (u> -> 7r°7)/S cxp (o; ->■ 37r) = 
(8.97 ± 0.16) x 10~ 2 [IT}, which is precisely the ratio we 
argued in Sec. IHII to be a pure prediction due to the 7T7r 
P-wave phase shift, independent of any subtraction con- 
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stant; for this quantity, our numerical result amounts to 
B(w -> 7r° 7 ) 



B(w 3tt) 



(8.39... 8. 69) x 10" 



(20) 



hence suggesting a saturation of the sum rule even at the 
95% level. 

We stress, however, that due to the slow convergence 
behavior of the integrand in Eq. (|12|l , we do not consider 
the sum-rule results to be extremely reliable: they de- 
pend rather strongly on the assumed intermediate and 
high-energy behavior of the tttt phase shift. For exam- 
ple, using a cutoff of A = 1.8 GeV in the dispersion inte- 
gral (|12p beyond which the asymptotic fall-off is enforced 
by hand, we find that this asymptotic region s > A 2 still 
yields a 10% correction to the uj — > 7r°7 branching ra- 
tio. We therefore rather take these as benchmark values 
to test the accuracy of the approximation of using only 
two-pion intermediate states in the dispersion relation: 
we expect this to work better in the description of the 
s-dependence of the transition form factor, in which we 
choose the subtraction constant in Eq. (TTTj) fixed to the 
experimental values of the V — > tt°j partial widths. The 
errors on these values contribute a large part to the un- 
certainty of the transition form factor and the differential 
V — > tt°£ + £~ decay width, which we will present in the 
following. 



A. w ->■ tt°£ + £- 

We start off by considering V — uj. In Fig. [5] we 
compare the absolute value squared of the uj — > 7r°7* 
transition form factor (left panel) and the differential 
uj —> tt°£ + £~ decay widths (right panel) calculated in our 
approach, standard VMD with a finite energy-dependent 
width, 3 and a chiral Lagrangian treatment with light vec- 
tor mesons from Refs. |34l. l35j. to data from Refs. 043 ■ 
The dispersive approach leads to a significant enhance- 
ment of the transition form factor over the pure VMD 
result, which in turn results in an improved description 
of the data. Part of this enhancement is even present if 
we use a simplified, VMD-inspired uj — ¥ 37r partial wave 
fi(s) = aft(s) inside the dispersion integral. As Fig. [3] 
suggests, using the correct full uj — > 37r P wave mainly 
leads to a further enhancement for invariant masses of the 
lepton pair near and slightly above the two-pion thresh- 
old. We note that using the slightly smaller sum-rule 
value for the normalization of the form factor (instead 
of the one determined from the experimental uj — > 7r°7 
width) would further enhance \F un o(s) — 1| 2 by 5-10%, 
albeit at the expense of a significantly enlarged uncer- 
tainty. However, we also find that our analysis cannot 
account for the steep rise towards the end of the decay 



region, which is somewhat better described by the calcu- 
lation in Refs. [H], HBj]. The size of the discrepancy for 
large invariant masses is surprising (note that the form 
factor in Fig. [5] is shown on a logarithmic scale) , in par- 
ticular given the level of agreement found in the sum 
rule for the uj — > 7r°7 branching ratio that should con- 
verge rather worse. Within the dispersive framework it 
is therefore hard to think of a plausible explanation for 
such a steep rise. We note that in contrast to cf> — > 3n, we 
have not yet been able to test the uj — > 3ir decay ampli- 
tude against experimental precision studies of the Dalitz 
plot, so a remaining deficit in our input for /i(s) cannot 
rigorously be excluded. Still, given the analogy to the 
<f> —> 3ir study in Ref . [l4| , it is implausible that this can 
account for the size of the difference. 

The transition form factors are often characterized by 
their slope at s — 0: 



dF v ^o (s) 



(21) 



We quote this slope in units of M p 2 (where we use M p = 
775.5 MeV), such that VMD suggests b^o = 1M~ 2 . 
Our dispersive analysis yields 



b un o = (1.41... 1.45) M- 



(22) 



therefore a significant enhancement with respect to the 
VMD value, yet not as large as the theoretical value 
found in Refs. [H], IHj], b^o 2M~ 2 , and significantly 
smaller than the experimental determinations 6 W7r o = 



(3.72 ± 0.10 ± 0.03) M; 



]■ K 



(3.73 ± 0.04 ± 



0.05) M~ 2 [9j. We note, however, that the latter experi- 
mental extractions are in principle model-dependent, as 
they rely on a monopole parameterization. 

In order to improve on the comparison of our form- 
factor description to the data, one may think of sub- 
tracting Eq. dTTJ) once more and treating the additional 
subtraction constant as a free parameter, at the expense 
of spoiling the high-energy behavior of the transition 
form factor. The difference between the once- and twice- 
subtracted representation amounts to an additive term 
Ab^o x s, and it is rather obvious that this term can- 
not account for the strong curvature in the form factor 
at higher energies, such that the overall picture is not 
drastically improved. Furthermore, the value given in 
Eq. (|2"21 amounts to a value for the slope given by a sum 
rule, which would be expected to converge much better 
than the one for / W7r o(0) in Eq. (fi"2"j) . yet it yields a result 
ostensibly off by a large factor. 

The differential decay width for uj — > 7r°/i + /i~ (top, 
right panel in Fig. [5]) is calculated according to Eq. 



3 This produces an almost negligible effect for V = ui, but guaran- 
tees sensible results for V = <fi. 



4 The normalization of the VMD prediction is obtained from the 
experimental uj — ¥ 7r°7 partial width, similar to the dispersive 
calculation. We refrain from displaying errors on the VMD cal- 
culation thus induced, since it merely serves illustrative purposes. 
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FIG. 5: Left panel: Numerical results for the <j — > 7r°7* transition form factor. Top, right panel: Differential lj — s- ir°n + fj,~ 
decay width. Bottom, right panel: Differential lu —¥ Tr°e + e~ decay width. Data for the transition form factor and the differential 
oj — s- TT°fi + fi~ width is taken from Refs. @-@] (we have not included the data set from Ref. [36l ] due to its fairly low statistics). 
We show pure VMD (dotted line), the results of a chiral Lagrangian treatment with explicit vector mesons [341 ] (white shaded 
curve with solid borders), and the dispersive solution for fi(s) = afi(s) (gray shaded curve with dashed borders) as well as the 
full dispersive solution (black hatched curve with solid borders). For lj — > Tv°e + e~ we do not display the pure Omnes solution, 
since it is virtually indistinguishable from the full dispersive result due to the strongly dominating kinematical factor in Eq. ((4} . 
The inset magnifies the region above the two-muon threshold (vertical dashed line). 



We observe that the values of the form factor close to the 
end of the decay region are actually strongly suppressed 
by phase space. From that vantage point the situation 
does not look as dire as when the form factor is considered 
directly; however, due to the smallness of the errors of 
those values our solution still deviates by several a. The 
integration of the spectrum yields 

B(u -> ttV + /-0 = (0.94 . . . 1.00) x 10~ 4 , (23) 

which agrees with the experimental average B cxp (uj —> 
ttV+aO = (1.3 ± 0.4) x 10~ 4 |3 within errors. This 
is not surprising: as the largest deviations from the ex- 
perimental form factor are strongly suppressed by phase 
space, they do not have a large influence on the partial 
width. 

We also display the ui — > 7r°e + e~ differential decay 
width (bottom, right panel in Fig. [5]), which has not been 
measured yet. Phase space combined with the 1 / s behav- 
ior of the virtual photon lead to a strong enhancement 
near threshold and a variation of the spectrum over many 



orders of magnitude; we therefore only display the full 
dispersive result, since it is almost indistinguishable from 
/i(s) = af2(s) on this scale. For a better comparison to 
lu — > 7T°/i + /!~, we also show this spectrum restricted to 
energies t/s > 2m^. As can be seen in Fig. [5j both lep- 
tonic final states yield very similar amounts of events in 
this energy range, where form-factor effects (deviations 
from pure QED) are felt most strongly. The integrated 
spectrum for lu — > 7r°e + e~ yields 



B(uj ->■ 7r°e + e-) = (7.6 . . . 8.1) x 10~ 4 , (24) 



where the uncertainty is dominated by the normalization 
given by B(ui — > tt°j) — the s-dependent e + e~ spectrum 
is largely given by pure QED. Equation ([24]) is in per- 
fect agreement with the experimental value B ckp (lu —> 
7r°e + e~) = (7.7 ± 0.6) x 10 -4 within uncertainties. 
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FIG. 6: Left panel: Numerical results for the (j> — ► it "f transition form factor. Top, right panel: Differential <j) — > ir n + /j,~ 
decay width. Bottom, right panel: Differential 4> — > ir°e + e~ decay width. We show pure VMD (dotted line), the dispersive 
solution for /i(s) = afl(s) (gray shaded curve with dashed borders), and the full dispersive solution with one subtraction 
(black hatched curve with solid borders) and two subtractions (white shaded curve with dashed-dotted borders) in the <j> —¥ 3n 
partial wave. For (f> — » Tv°e + e~ we only show the once-subtracted dispersive solution, since again neither the Omnes solution 
nor the twice-subtracted one is visibly distinguishable from the once-subtracted result on the scale shown. The inset magnifies 
the region above the two-muon threshold (vertical dashed line). 



B. tj>-t tt°1 + £- 

Next we consider V — 4>. The results are displayed in 
Fig. [B] for the absolute value squared of the cj) —> 7r°7* 
transition form factor (left panel) and the differential 
4> — > ir°£ + £~ decay widths (right panel). There are no 
experimental data in any of the decay channels yet. For 
the <f> — > 7r°7* transition form factor we observe a sim- 
ilar behavior as for the cj: in the full decay region, the 
form factor calculated with the dispersive approach is en- 
hanced compared to the pure VMD result; in addition, 
we observe the two-pion-threshold enhancement of the 
full dispersive result with respect to fi(s) = afi(s). Due 
to the strong rise of the full solution for /i (s) towards this 
threshold, see Fig. [3J the corresponding F^o (s) almost 
approaches a cusp-like behavior at s = . Since the <f> 
as opposed to the cj transition form factor encompasses 
the p resonance region, we can also observe that the full 
solution for /i(s) slightly reduces the height of the res- 
onance peak with respect to the simplified assumption 
fl(s) = aO(s), which agrees with our observations in 



Ref. |14[. We note that using a twice-subtracted disper- 
sion relation for the partial- wave amplitude fi(s), with 
the additional subtraction constant fitted to the <p — > 37r 
Dalitz plot of Ref. [32[, does not change our results by 
all that much: the differences are smaller than the over- 
all uncertainty in our transition-form-factor prediction. 
This corroborates our skepticism that an imperfect de- 
termination of the cj — v 37r P wave is the likely source 
of the discrepancy seen in the cj transition-form-factor 
data. 

Again, we also quote the derivative of the form factor 
at s = 0: 

b^o = (1.52... 1.61) M- 2 , (25) 

which is still somewhat larger than b U7r o, see Eq. (|22[) . 
but again not nearly as large as the slopes found experi- 
mentally in other vector-meson conversion decays. 

The observations above concerning the differences of 
the various theoretical predictions translate directly to 
the 4> — ¥ 7r°/z + /i _ differential decay spectrum (top, right 
panel of Fig. [5]). We find that the p resonance leaves a 
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clear imprint on the spectrum, as one observes a second 
peak structure that counterbalances the drop-off of the 
phase-space factor. The integrated spectrum yields 

B oncc ((f> -> 7T V/O = (3-7 . . . 4.0) x 1(T 6 , 
gtwice^ ^ tt V/O = (3.8 . . . 4.1) x 1CT 6 , (26) 

for the once- and twice-subtracted 4> — > 37r partial-wave 
amplitudes, respectively, perfectly compatible within the 
error ranges. There is currently no experimental mea- 
surement of the partial width to be compared with. 

As for the corresponding o> decay, the differential (f> — ¥ 
ir°e + e~ decay width is enhanced for small s by several or- 
ders of magnitude; for this reason, we only display the full 
dispersive solution based on the once-subtracted (f> 3tt 
partial wave /i(s), the alternatives being indistinguish- 
able on this scale. Again, an insert concentrates on ener- 
gies above the two-muon threshold for better comparison 
of the expected event rates in both final states. The re- 
sults for the integrated spectra are 

gonce^ ^ ^ e + e -) = (1.39 . . . 1.51) x 1Q-5 , 
gtwice^ ^Og+g-) = (i 40 _ _ _ 153) x i -5 ; ( 27 ) 

for both of the full solutions, respectively. Compared 
with the experimental value of B c * p ((f> —> ir°e + e~) = 
(1.12 ± 0.28) x 10~ 5 [H, we find agreement within un- 
certainties. 

We wish to emphasize the significance of an experi- 
mental investigation of the <f> — > 7r°7* transition form 
factor. Deviations from the VMD picture now seem to 
be well established in lu — > 7r°7*; strikingly enough, both 
this and the latest measurement of the transition form 
factor in (f> — > 777* (Tlj . when parameterized in terms 
of a monopole form factor, yield monopole mass param- 
eters significantly below the scale of the physical vec- 
tor mesons, but (of course) too large to be accessible 
within the physical decay region. This is different in 
the p resonance can be measured in this de- 
cay; if there systematically is a steep form factor rise as 
seen in u> — > 7r°7*, mapping it out in full in 4> — > 7r°7* 
will help clarify its origin. From the theoretical side, our 
dispersive analysis for this process is based on a very 
precisely measured <^> — > 37r Dalitz plot, such that we 
are very confident about the reliability of our prediction. 
We thus strongly advocate an experimental analysis of 
the 4> — > 7r°7* form factor to the best possible precision. 

As a final illustration, we show the (experimentally 
unobservable) phases of the V — > 7r°7* transition form 
factors, both for w and (f>, in Fig. [7J We calculate these 
from the unsubtracted solution to the discontinuity equa- 
tion (|10l) , as we can only fix the modulus of the subtrac- 
tion constant fy^o (0) by means of the V — > 7r°7 partial 
width, not its phase, which is nonvanishing due to the 
complex discontinuity of the V — > 37r partial- wave ampli- 
tude fi(s). Only the phases of the full dispersive calcula- 
tions, as compared to the irir P-wave phase shift, are dis- 
played. We refrain from showing the (small) error bands 




[GeV] 



FIG. 7: Phase of the electromagnetic transition form factor 
for cj — > 7r°7* (solid line) and cj> — ► 7r°7* (dotted line) to- 
gether with the input phase-shift parameterization (dashed 
line). The inset shows a magnification of the region below 
the two-pion threshold. 




FIG. 8: Three-pion-cut contribution to the V — > tv°£ + £ tran- 
sition vector form factor, inducing an imaginary part also for 
s < AMI. 

and fix the input in complete analogy to Fig. [3l An ad- 
ditional consequence of the complex discontinuity of the 
V ^ 3tt partial-wave amplitude is that Watson's final- 
state theorem also does not apply to the transition form 
factors, and their phases are different from 8{s): three- 
pion-cut effects, see Fig. [8j produce nonvanishing imag- 
inary parts/nonvanishing phases of the transition form 
factors also below the 7T7t threshold, s < 4M%. We ob- 
serve the transition- form-factor phases to be significantly 
larger than S(s) above the two-pion threshold, and a ten- 
dency to small negative values below. 



V. SUMMARY 

In this article, we have analyzed the V — > 7r°7* elec- 
tromagnetic transition form factors for V — to and V = <f> 
by means of a dispersive framework. It requires the cor- 
responding P-wave projection of the V — > ir + TT~ir° decay 
amplitudes, and the pion vector form factor as input, 
both of which depend on the pion-pion P-wave scattering 
phase shift as input and are otherwise predictions up to a 
subtraction constant determining the overall normaliza- 
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tion of the V — > 3ir amplitudes. The consistent treatment 
of crossed-channel effects in the V — > 3ir partial-wave 
amplitudes by incorporating three-particle cuts leads to 
a nontrivial analytic structure for the transition form fac- 
tors, in particular its phase does not follow the tttt P-wave 
phase. 

We have calculated the real-photon V — > tt°j branch- 
ing ratios using a sum rule, which yields good agreement 
with the experimental 4> ~^ 7r °7 branching ratio and in- 
dicates that the sum rule for uj — > 7r°7 (which is much 
more precisely determined experimentally) is saturated 
roughly at the 90%-95% level by two-pion intermediate 
states. To lessen the dependence on medium-to-high- 
energy input, we have oversubtracted the dispersion re- 
lation for the form factors and used the real-photon par- 
tial widths as input for the subtraction constant. We 
found that this approach leads to an enhancement com- 
pared to a pure VMD calculation and thus to an im- 
proved description of experimental data from NA60 for 
the lu — > 7T°/i + /!~ channel. Three-pion effects in partic- 
ular lead to an enhancement in the two-pion-threshold 
region. 

We are unable to solve the puzzle of the steep rise in 
the w — > 7r°7* transition form factor data close to the 
end of the decay region. In order to try to better under- 
stand the physical mechanism behind this enhancement, 
we strongly advocate a measurement of the <fi ~~ * 7r°7* 
transition form factor: the fact that the physical region 
of the decay goes beyond that of the corresponding u de- 
cay and incorporates the p resonance peak suggests that 
it should give some clues about the nature of this rise. 

While our predictions for branching ratios of the var- 
ious V —> tt°£ + £~ channels are in good agreement with 
experimental determinations, data on decay spectra only 
exists for u — » Tr°p + p~. It would certainly be helpful, 
especially in light of a theoretical analysis of contribu- 
tions to light-by-light scattering, if precision data for ad- 
ditional channels could be obtained [37] . 
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Appendix A: Pion vector form factor representation 
including higher resonances 

To account for the effects of higher resonances in the 
pion form factor, at least in the elastic approximation, 
we use the analytic representation [29j 



* [S) Ml - s - iM p T p {s) 



exp 



7 s e 



M%,-s-iM p ,T p ,(s) 



exp 



exp 



96tt 2 F2 
sr p /A 7r (g) 

ttM^M 2 ,,) j ' 
(Al) 



where 



Ms) = log ^ + ~f ~ 3 + (a) log , 



<(s), <T„(S) = 4/1 



P'.P" 



96ttF2 
M n , n „ ( s - AMI 



.3/2 



p',p" i 



(A2) 



and is the pion decay constant. We fit the various 
masses, widths, and coupling constants to the t~ — > 
■k^h°v t data of Ref. [26[ (thus eschewing the additional 
complication of p—ui mixing). Note that we have omit- 
ted the kaon loop contributions to T p (s), so that our fit 
values are slightly different as compared to Ref. [29| . We 
obtain 



We are grateful to S. Leupold for extensive, very use- 
ful discussions and a critical reading of this manuscript. 
We would like to thank S. Damjanovic and G. Usai for 
providing us with the data of Refs. [1, and G. Colan- 
gelo for the preliminary 7T7T phase-shift results of Ref. [24| . 
Partial financial support by the DFG (SFB/TR 16, "Sub- 



Mp> = 1.44 ±0.01 GeV 



7 = 0.097 ±0.009 



M n 



0.34 ±0.03 GeV , 
0.5 ±0.2 , 



1.71 ±0.05 GeV , T p » — 0.13 ± 0.03 GeV , 



-0.02 ±0.02 , 



1.1 ±0.6 



(A3) 
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